Hawking radiation from dynamical horizons 
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In completely local settings, we establish that a dynamically evolving black hole horizon can be 
assigned a Hawking temperature. Moreover, we calculate the Hawking flux and show that the radius 
of the horizon shrinks. 
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PACS numbers: 

The laws of black hole mechanics in general relativity 
are remarkably analogous to the laws of thermodynam- 
ics This analogy is exact when quantum effects are 
taken into account. Indeed, Hawking's semiclassical anal- 
ysis establishes that quantum mechanically, a stationary 
black hole with surface gravity k radiates particles to in- 
finity with a perfect black body spectrum at temperature 
k/2-k 0. Consequently, asymptotic observers perceive a 
thermal state and assign a physical temperature to the 
black hole. The precise match to thermodynamics is com- 
plete when the thermodynamic entropy of the black hole 
is identified with a quarter of its area [3[ . 

The original calculation of Hawking is independent of 
the gravitational field equations. It relies only on the 
behavior of quantum fields in a specific spacetime geom- 
etry describing a stationary black hole formed due to a 
gravitational collapse. Over the years, several other tech- 
niques have been developed to study spontaneous particle 
emission and the Hawking temperature for more general 
spacetimes. For example, the Hartlc-Hawking proposal 
@ and the Euclidean approach @ have been extensively 
used to associate thermal states to spacetimes with bi- 
furcate Killing horizons. In fact, it has been established 
that in any globally hyperbolic spacetime with bifurcate 
Killing horizon, there can exist a vacuum thermal state 
at temperature k/2tt which remains invariant under the 
isometries generating the horizon Q . 

Although these constructions are elegant, they arc 
quite restrictive, inapplicable even for spacetimes with 
superradiance [7[. These formulations also do not indi- 
cate how such a thermal state may arise as a result of 
some version of physical process. In addition, their ex- 
istence requires knowledge of global structure of space- 
time. As a result, they do not appear very useful to study 
thermal properties of local horizons. On the other hand, 
the laws of black hole mechanics apply equally well to 
black hole horizons which can been proved using only lo- 
cal geometrical properties of null surfaces, without any 
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assumptions on the global development of the spacetime 
in which the horizon is embedded [cj-fllT] . It has also been 
established that such horizons can be assigned an ent ropy 
proportional to the area of the local horizon [12, Il3j |. 
Thus, it seems to be a reasonable physical expectation 
that even with a local definition of black hole horizon 
one should be able to establish the analogy to thermo- 
dynamics. More precisely, such horizons should have a 
temperature of k/2tt. 

Incidentally, this question has been investigated in a 
semiclassical approach which treats Hawking radiation as 
a quantum tunneling phenomenon [Hj]-[l8j]. The method 
involves calculating the imaginary part of the action for 
the (classically forbidden) process of s-wave emission, 
from inside and through the horizon (see fl9j for more 
details). Using the WKB-approximation the tunnelling 
probability for such a classically forbidden trajectory is 
calculated to be, V = e ~ 21mS where, S is the classical 
action of the trajectory to leading order in h. This is 
equated to the Boltzmann factor e~^ E to extract the in- 
verse Hawking temperature /3. The main advantage of 
this formalism is that the calculations involve only the 
local geometry and hence can be applied to any local 
horizon. Indeed, tunnelling method has been applied to 
local dynamical black hole horizons and the temperature 
is found to be k/2tt where k is the dynamical surface 
gravity [2(| HH . Still there are some problems with the 
method itself and some issues which have not been ad- 
dressed in this treatment of dynamical horizons. First, 
the approach depends heavily on the semiclassical ap- 
proximation and though it is argued that this remains 
valid near the horizon, it would be better to devise a 
more general formalism which does not rely on WKB-likc 
approximations. Secondly, in calculating the imaginary 
part of the semiclassical action S from the Hamilton- 
Jacobi equation, a singular integral appears with a pole 
at the horizon. While for the static case the result is stan- 
dard, for the dynamical horizon it is not clear how the 
integration is to be performed since the position of the 
horizon changes in a dynamical process. Lastly, in all 
these treatments of radiation from dynamical horizons 
the evolution of the horizon itself is never addressed. In 
other words, it is not clear how the horizon loses area due 
to emission of a flux of radiation. The local formalisms of 
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black hole horizon should be able to address these issues. 

In this paper, a formalism is developed to establish two 
basic issues. First, that one can associate a temperature 
to local dynamical horizons without the need for any 
WKB-like approximation schemes. Second, that there 
exists a precise relation between the radiation emitted 
by the horizon and area loss, i.e., flux of outgoing radia- 
tion through the horizon in between two partial Cauchy 
slices exactly equals the difference of radii of the sphere 
that foliates the horizon at those two instances. 

We elucidate our arguments as follows. To calculate 
temperature for local dynamical horizons, we begin by 
considering the Kodama vector field [22[ • For dynamical 
spacetimes, this vector field provides a preferred timelike 
direction and is parallel to the Killing vector at spatial 
infinity which we assume to be fiat. We can construct 
well-behaved positive frequency field modes on both sides 
of the horizon by considering the Kodama vector field 
but the outgoing modes exhibit logarithmic singularities 
on the horizon under some approximation. However, if 
considered as distribution valued, these modes can be 
interpreted as horizon crossing and the probability cur- 
rent for these modes remain well defined. The Hawking 
temperature is determined if one equates the conditional 
probability, that modes incident on one side is emitted 
to the other side, to the Boltzmann factor [23l . |24| , 
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Since this method does not depend on the entire evo- 
lution of the field modes in the spacetime, it is ideally 
suited for our purpose. 

To evaluate the Hawking flux, we recall that there 
are two well known (and related) definitions of local 
black hole horizon, the future outward trapping horizon 
(FOTH) HH and the Dynamical Horizon [2! [2^ (or its 
equilibrium version called the isolated horizon) . In these 
local settings, black hole horizons are a stack of appar- 
ent horizons which, under suitable energy conditions, are 
either null or spacelike. As such, energy flux can only 
remain on the surface or flow into the horizon. In or- 
der that matter fields flow out off such a surface requires 
that the surface must be timelike in some afflne interval. 
However, to achieve a timelike evolution of the horizon, 
some energy conditions need to be violated. This is only 
natural since Hawking radiation necessarily associates, 
with the thermal emission of particles, a positive flux of 
energy flowing to infinity (we shall assume that the space- 
time is asymptotically flat) and a corresponding flux of 
negative energy flowing into the black hole (this negative 
energy flux can also be motivated by the fact that the ex- 
pectation values of stress energy tensor of quantum fields 
generically violate energy conditions). In this process the 
horizon looses area and energy. 

The plan of the paper is as follows: First, we will dis- 
cuss the geometrical setup which is based on future outer 
trapping horizon (FOTH). Next, we show that how the 
Hawking temperature is proportional to the dynamical 



surface gravity associated with the Kodama vector. Fi- 
nally, we will calculate the flux of energy radiated in a 
dynamical process. 

We begin with definitions. We follow the conventions 
of @. Consider a four dimensional spacetime ^# with 
signature (—,+,+,+). A three-dimensional submanifold 
A in ^# is said to be a future outer trapping horizon 
(FOTH) if 1 ) It is foliated by a preferred family of topo- 
logical two-spheres such that, on each leaf S, the expan- 
sion 6 + of a null normal vanishes and the expansion 
0_ of the other null normal l°_ is negative definite, 2) 
The directional derivative of 9 + along the null normal 
(i.e., «Sfj _ 6* + ) is negative definite. 

Thus, A is foliated by marginally trapped two-spheres. 
According to a theorem due to Hawking, the topology of 
S is necessarily spherical in order that matter or gravita- 
tional flux across A is non-zero. If these fluxes are iden- 
tically zero then A becomes a Killing or isolated horizon. 

Even though our arguments will remain local, for def- 
initeness, we choose a spherically symmetric background 
metric 

ds 2 = -2e~ f dx + dx~ + r 2 (d8 2 + sin 2 6dcj) 2 ) (2) 

where both / and r are smooth functions of x ± . The ex- 
pansions of the two null normals arc 6± = (2/r) d±r re- 
spectively where d± = d/dx ± . In this coordinate system, 
the second requirement for FOTH translates to d-0 + < 
on A. 

Let the vector field t a = 1°: + hit be tangential to the 
FOTH for some smooth function h. Then the Raychaud- 
huri equation for Z5. and the Einstein equation implies 

d+9+ = -hd-9 + = -8itT ++ . (3) 

where T ++ = T a b l+l b + and T a t is the energy momentum 
tensor. Several consequences follow from this equation. 
First, the FOTH is degenerate (or null) if and only if 
T ++ = on A. In that case, the FOTH is generated 
by C. Degenerate FOTH is not interesting for Hawk- 
ing radiation because this implies d+r = 0. As a conse- 
quence, the area, A = Airr 2 of S, and the Misner-Sharp 
energy for this spacetime, given by E = |r, also remains 
unchanged. Secondly, since t 2 = —2he~f, a FOTH be- 
comes spacclikc if and only if T++ > and is timclikc if 
and only if T ++ < 0. 

For a timclikc FOTH, several consequences follow. 
Here, „S? t r < 0, and hence, A is timclikc if and only if the 
area A and the Misner-Sharp energy E decreases along 
the horizon. This is also expected on general grounds 
since the horizon receives an incoming flux of negative 
energy, T ++ < 0. 

As we have emphasized before, in the dynamical space- 
time the Kodama vector field plays the analog role of 
the Killing vector. For this spacetime, it is given by 

K a = e f (d-r) d% - e f (d+r) d a _ . (4) 

The surface gravity is defined through AT a V^A' a j = kK), 
and is k = -e f d-d+r. The FOTH condition d-6+ < 
implies k > 0. 
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Let us now determine the positive frequency modes 
of the Kodama vector. It is easy to see that any smooth 
function of r is a zero-mode of the Kodama vector. Once, 
a zero-mode is obtained, other positive frequency eigen- 
modes are evaluated using 



iK Z^= ujZ^ 



(5) 



Here, Z^ are the eigenfunctions corresponding to the pos- 
itive frequency u>. For simplification, let us introduce 
new coordinates, y = x~ and r and two new functions, 
Z u (y,r) = Z UJ (x + 1 x^) and G(y,r) = e f (d + r). As a 
result, the eigenvalue equation reduces to 



GdyZu = iuj Z^ 



(0) 



Integrating and transforming back to old coordinates, the 
above equation gives 



= F{r) cxp (iuj J 



dx 



■fd 4 



(7) 



where F(r) is an arbitrary smooth function of r and the 
subscript r under the integral sign denotes that while 
doing the integration r is kept fixed. To evaluate the 
integral in (|7|), we multiply the numerator and the de- 
nominator by {8-6+) and use the fact that for any fixed 
r surface, (<9_0+) = — 2k/r, (although the strict in- 
terpretation of k as the surface gravity holds only for 
surfaces with 9 + = 0, it exists as a function in any neigh- 
bourhood of the horizon). Thus, in some neighbourhood 
of the horizon we get 



dx d_Q + 
ef d+rd-9. 



d9+ 
k9~T 



(8) 



We now assume (this is the only assumption we make in 
this calculation) that during the dynamical evolution k is 
a slowly varying function in some small neighbourhood of 
the horizon (the zeroth law takes care of it on the horizon, 
but we also assume it to hold in a small neighbourhood 
of the horizon). This gives 



Zuj — 




»+ < o. 



(9) 



where the spheres are not trapped 'outside the trapping 
horizon' (9+ > 0) and fully trapped 'inside' (9+ < 0). 
These are precisely the modes which are defined outside 
and inside the dynamical horizon respectively but not on 
the horizon. Now we have to keep in mind the modes ((9]) 
are not ordinary functions, but are distribution-valued. 
Comparing with the spherically symmetric static case 
[H, we find for e -> 0+ 



for 9+ > 
for 9+ < 



(10) 



for the choice A = —iui/k. The distribution (|TU1) is well- 
defined for all values of 9+ and A, and it is differentiablc 



to all orders. The modes Z* are given by the complex 
conjugate distribution. 

We wish to calculate the probability density in a sin- 
gle particle Hilbcrt space for positive frequency solutions 
across the dynamical horizon 



Z,.,KZu) — KZ,.,Zu 



= ioZ*Z u 



[11) 



A straightforward calculation gives, apart from a positive 
function of r, 



g(w) = uj(9 + + ie) k (9+ - ie)~^ 

juj for 9+ > 

I ueT for 9 + < 0. 



(12) 



The conditional probability that a particle emits when it 
is incident on the horizon from inside is, 



cm,isston\ inci 



dent' 



= P k 



(13) 



This gives the correct Boltzmann weight with the tem- 
perature fc/271", which is the desired value. 

We now show that as the horizon evolves, the ra- 
dius of the 2-sphere foliating the horizon shrinks in pre- 
cise accordance with the amount of flux radiated by the 
horizon. To study the flux equation, consider new co- 
ordinates, (x + ,x~) i Y (9 +1 x~) where x~ = x~ . On 
FOTH, (<9_0 + )/(<9 + 0+) is equal to -(<9_<9 + r)/(47rr T ++ ) 
and negative definite. As a result, the derivatives are 
related to each other by 



( \ d + . 



Amr T. 



++ 



(14) 



It is not difficult to show that cL is proportional to the 
tangent vector t a to the FOTH. Observe that the normal 
one- form to A must be proportional to (dr — 2 dE). which 
on the horizon is equal to the one-form 



lire f r 2 T ++ - 2re f d_d + r) <9_r dx' 



(15) 



In arriving at the above identity we have made use of two 
Einstein's equations @ 

rd-d+r + d+r d^r + \e~ f = 47rr 2 T_ + , (16) 
d 2 + r + d+Jd+r = -A-KrT ++1 

and energy equations 

d±E = 27re / r 3 (T_ + ± -T ±± fl T ). (17) 

As a result, the normal vector n a is proportional to 



d+ 



Aitr T. 



++ 



<9_<9 4 



d- = d+ — hd- 



(18) 



so that the tangent vector t a = di_ + hd°L , which is clearly 
proportional to (fT4")l . 
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So x ,0,cj) are natural coordinates on FOTH. The line- 
element © induces a line-element on A 



ds 2 = -2e- / /i- 1 (d.i-) 2 + r 2 (d9 2 + 



sm 



; ). (19) 



Consequently, the volume element on A is given by 
dfx = \ / 2e~fh~ 1 r 2 sinO dx~ d6d<f>. We can now calcu- 
late the flux of matter energy that crosses the dynamical 
horizon — it is an integral on a slice of horizon bounded 
by two spherical sections S\ and S2 



■9 



dfi T ab h a K b 



where n a is the unit normal vector 
1 



V2he-f 



f {d a + -hd a _) 



(20) 



(21) 



and K a is the Kodama vector. Using spherical symmetry, 
cqn. (fT4|) and eqn. (fl"6)) . we get 



& = J dx~ 4:TTr 2 (^-T ++ - T + -)e f d^r 

= I dx~ 47rr 2 (-^-<9 + c^r-T + _)e / <9_r. (22) 
J iirr 

Making use of the Einstein equation (fTS)) on the horizon 
and (fl"4|) . we get 



dx —d-T 



:(r-2 - n) 



dx 



(23) 



where r± , ri are respectively the two radii of Si, 1S2. Since 
the area is decreasing along the horizon, r<i < r± where 
S2 lies in the future of Si. As a result, the outgoing flux 
of matter energy radiated by the dynamical horizon is 
positive definite (and the ingoing flux of matter energy 
is negative definite). The flux formula (f2"3")l differs from 
that given in [2|| . Since the Kodama vector field provides 
a timelike direction and is null on the horizon, it seems 
more appropriate to use K a for the dynamical horizon. 



The derivation of Hawking temperature and the flux 
law depends on two assumptions. First, that the Ko- 
dama vector exists in the spacetime. For spherically sym- 
metric spacetimes, the Kodama vector field exists unam- 
biguously and the Misncr-Sharp energy is well defined. 
For more general spacetimes, a Kodama-like vector field 
is not known, however, one can still define some mass 
for such cases that reduces to the Misncr-Sharp energy 
in the spherical limit (27j . The second assumption, the 
existence of a slowly varying k can also be motivated 
for large black holes. In such cases, the horizon evolves 
slowly enough so that the surface gravity function should 
vary slowly in some small neighbourhood of the horizon. 
Alternatively, we can conclude that the Hawking temper- 
ature for a dynamically evolving large black hole is fc/27r 
if the dynamical surface gravity is slowly varying in the 
vicinity of the horizon. 

The set-up described in this paper can be further devel- 
oped to model dynamically evaporating black hole hori- 
zons through Hawking radiation, analytically as well as 
numerically. Over the years, several models have been 
constructed which study radiating black holes, formed in 
a gravitational collapse, based on the imploding Vaidya 
metric with a negative energy-momentum tensor, show 
that a timelike apparent horizon forms due to violation of 
energy conditions [28j]. However, such models are based 
on global considerations of event horizons, while local 
structures like that used in [29j might be useful for a 
better understanding of Hawking radiation and compu- 
tations of quantum field theoretic effects (see also [30] ) - 

It is also interesting to speculate on the extension of 
the present method for other diffcomorphism invariant 
theories of gravity. While the zeroth and the first law 
hold for any arbitrary such theory, the second law has 
only been proved for a class of such theories [3l[. If the 
present formalism can be extended to other theories of 
gravity, it will lend a support to the existence of the area 
increase theorem for such theories. 

While more interesting and deeper issues can only 
be understood in a full quantum theory of gravity, the 
present framework can elucidate the suggestions of [32j 
and provide a better understanding of the Hawking radi- 
ation process. 



[1] J. M. Bardeen, B. Carter and S. W. Hawking, Commun. 

Math. Phys. 31, 161 (1973). 
[2] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975) 
[3] J. D. Bekenstein, Phys. Rev. D 7, 2333 (1973). 
[4] R. M. Wald, Chicago, USA: Univ. Pr. (1994) 205 p 
[5] J. B. Hartle and S. W. Hawking, Phys. Rev. D 13, 2188 

(1976). 

[6] G. W. Gibbons and M. J. Perry, Proc. Roy. Soc. Lond. 

A 358, 467 (1978). 
[7] B. S. Kay and R. M. Wald, Phys. Rept. 207, 49 (1991). 
[8] S. A. Hayward, Phys. Rev. D 49, 6467 (1994). 
[9] S. A. Hayward, Class. Quant. Grav. 15, 3147 (1998) 



[arXiv:gr-qc/9710089] . 
[10] A. Ashtekar, C. Beetle, O. Dreyer, S. Fairhurst, B. Kr- 

ishnan, J. Lewandowski and J. Wisniewski, Phys. Rev. 

Lett. 85, 3564 (2000) [arXiv:gr-qc/0006006j . 
[11] I. Boot h and S. Fairhurst, P hys. Rev. Lett. 92, 011102 

(2004) [arXiv:gr-qc/0307087| . 
[12] A. Ashtekar, J. Baez, A. Corichi and K. Krasnov, Phys. 

Rev. Lett. 80, 904 (1998) [arXiv:gr-qc/9710007| . 
[13] A. Gho sh and A. Perez , Phys. Rev. Lett. 107, 241301 

(2011) |arXiv: 1107. 13201 [gr-qc]]. 
[14] M. K. Parikh and F. Wilczek, Phys. Rev. Lett. 85, 5042-5 

(2000). 



5 



[15] S. Shankaranarayanan, T. Padmanabhan and K. Srini- 
vasan, Class. Quantum Grav. 19, 2671-88 (2002). 

[16] M. Angheben, M. Nadalini, L. Vanzo and S. Zerbini, J. 
High Energy Phys. 0505 014 (2005). 

[17] R. Kerner and R. B. Mann, Phys. Rev. D 73 104010 
(2006). 

[18] B. Chatterjee, A. Ghosh and P. Mitra, Phys. Lett. B 661 
307-11 (2008). 

[19] L. Vanzo, G. Acquaviva and R. Di. Criscienzo Class. 

Quantum Grav. 28 183001 (2011). 
[20] R. Di. Criscienzo, M. Nadalini, L. Vanzo, S. Zerbini and 

G. Zoccatelli, Phys. Lett. B 657 107-11 (2007). 
[21] S. A. Hayward, R. Di Criscienzo, L. Vanzo, M. Nadalini 

and S. Zerbini, C lass. Quant. Grav. 26, 062001 (2009) 

|arXiv:0806.0"0l4l [gr-qc] ] . 
[22] H. Kodama, Prog. Theor. Phys. 63, 1217 (1980). 
[23] T. Damour and R. Ruffini, P hys. Rev. D 14, 3 32 (1976). 
[24] B. Chatterjee and A. Ghosh. larXiv:1201.4017l [gr-qc], 
[25] A. Ashtekar and B. Krishnan, Phys. Rev. Lett. 89, 



261101 (2002) [arXiv:gr-q c/0207080 ; Phys. Rev. D 68, 

104030 (2003) |arXiv:gr-qc/ 0308033 . 
[26] A. Ashtekar and B. Krishnan, Living Rev. Rel. 7, 10 

(2004) [arXiv:gr-qc/0407042| . 
[27] S. Mukohyama and S. A. Hayward, Class. Quant. Grav. 

17, 2153 (2000) [arXiv:gr-qc/9905085l . 
[28] W. A. Hiscock, Phys. Rev. D 23, 2813 (1981); R. Bal- 

binot, R. Bergamini and B. Giorgini, Nuovo Cim. B 70 

(1982) 201. 

[29] S. A. Haywar d, Phys. Rev. Let t. 96, 031103 (2006). 

[30] D. Pranzetti. larXiv: 1204.07021 [gr-qc], 

[31] A. Chatterjee and S. Sarkar, Phys. Rev. Lett. 108, 
091301 (2012) [arXivillll. 30211 [gr-qc]]; S. K olekar, 
T. Padmanabhan and S. Sarkar, arXiv:1201.2947 [gr-qc]. 

[32] A. Ashtekar and M. Bojowald, Class. Quant. Grav. 
22, 3349 (2005) [arXiv:gr-qc/0504029l ; A. Ashtekar, 
V. Taveras and M. Varadarajan, Phys. Rev. Lett. 100, 
211302 (2008) |arXiv:0801.18TT1 [gr-qc]]. 



